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The random walk of an electrostatic field using parallel infinite charged planes 
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We show that it is possible to generate a random walk with an electrostatic field by means of several parallel 
infinite charged planes in which the surface charge distribution could be either zbcr. We formulate the problem 
of this stochastic process by using a rate equation for the most probable value for the electrostatic field subject 
to the appropriate transition probabilities according to the electrostatic boundary conditions. Our model gives 
rise to a stochastic law when the charge distribution is not deterministic. The probability distribution of the 
electrostatic field intensity, the mean value of the electrostatic force and the energy density are obtained. 
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I. INTRODUCTION 

In a typical electrostatic problem you are given a volume 
charge distribution p(r) and you want to find the electric field 
E it produces. The fundamental relation between the source 
of the field, i.e. charge distribution, and the electric field is 
given by Gauss’s law in 

V-E=-^ (1) 

^0 

where cq is the electric permittivity of free space. For the par¬ 
ticular case of infinite charge sheets the electric field depends 
on one variable only, i.e. E(^), and the field strength it pro¬ 
duces is constant for all distances from the plane. The elec¬ 
tric field caused by an infinite sheet of volume charge density 
p{z) = (75{z) is given by O 

E = 7-sign(2:) (2) 

where a is the surface charge density and the function sign is 
defined as 16 ) 

. . . / +1 if 2 : > 0 

sign( 2 ) = |_^ 

and has been introduced to account for the vector nature of 
the electric field. The total electric field due to N of these 
infinite charge planes is then just the sum of the individual 
contributions, i.e. 

E = X (4) 

n=l 

where cr(zn) represents the surface charge density of the n-th 
charge plane at position Zn. Note that it is crucial to know 
the surface charge density on each infinite sheet to completely 
determine the value of the electric field through Eq. 0, if the 
source of the electric field is not known then we would need 


to approach the problem in a different way. We can formu¬ 
late this same problem from a probabilistic point of view us¬ 
ing random walks. The connection between electrostatics and 
Brownian motion or random walk was first shown in 1944 by 
Kakutani (T). Since then many investigators have applied the 
probabilistic potential-theory to solve differential and integral 
equations for calculating electrostatic field and capacitance in 
integrated circuits iliD. More recently, the connection be¬ 
tween electrostatics and quantum mechanics was first shown 
by one of the authors OlIID. The purpose of this paper is 
to give a description of the electrostatic field when the source 
of the field is not known. If the volume charge density is not 
known then we can not use Eq. Q to solve for the electric 
field, instead we need to use the boundary conditions for the 
electric field. For the restricted case of electrostatic charges 
an fields in vacuum, the appropriate boundary conditions for 
the electric field across a surface charge distribution is given 
by im 

^ahove '^helow ~ ^ ( 5 ) 

^0 

where n is a unit vector perpendicular to the surface. For the 
case of parallel infinite charged planes the boundary condition 
is given by 

^i^^above ^{^^below — ( 6 ) 

Equation ([^ says that, for the special case in which the infinite 
charged planes have a surface charge distribution given by a, 
the future electric field value depends only on the present elec¬ 
tric field value. This means that given precise information of 
the present state of the electric field value, the future electric 
field value does not depend on the past history of the process. 
This point is at the very heart of the Markov chain processes. 
This means that we can use a probabilistic model which gives 
a complete description of the electric field when the source of 
the field is not known. 

The problem that we would like to address in this paper is the 
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one where you have several parallel infinite charged planes 
placed along the z axis in which the surface charge density 
could be either ±cr. Then the field would evolve along the 
2: axis making random jumps each time it crosses an infinite 
charge sheet. We are going to show that the behavior of the 
electric field can be studied like a random walk. The motiva¬ 
tion of this paper is due to the fact that many physical systems 
display a non-deterministic disorder. A better understanding 
and prediction of the nature of these systems is achieved by 
considering the medium to be random. The randomness mod¬ 
els the effect of impurities on a physical system or fiuctua- 
tions. Randomness appears at two levels in our problem. It 
comes in the description of the intensity of the electric field 
and also comes in the description of the medium. 

The article is organized as follows. In section II we give 
a brief review of Markov chains. Then in section III we will 
apply the Markov chain theory to the case where we have sev¬ 
eral infinite charged planes placed along the z axis in which 
the surface charge density can be either ±cr. In section IV we 
present the numerical results and use our model to calculate 
the mean electrostatic force and energy density for our prob¬ 
lem. In the last section we summarize our conclusions. 


II. MARKOV CHAINS 

In this section we will give a brief explanation of the 
Markov chain theory concepts that we need for the formu¬ 
lation of the problem. To introduce these concepts lets imag¬ 
ine a random experiment like throwing a normal dice. In this 
case we have six possible outcomes which are 1, 2,3,4, 5 or 
6. We will denote the set of the possible outcomes with the 
letter S and name it as the state space. These outcomes de¬ 
fine a random variable X which can take values from S. Each 
outcome can be associated with their respective probability p 
which must take values only in the interval 0 < p < 1. We 
can write P{X = i) = Pi which means that the probability 
of the random variable X to take a value i its equal to pi with 

i G s'ini. 

Also if the same experiment is repeated n times, and in gen¬ 
eral each trial is dependent of the previous trials then Pi{n) 
represents the state of the random variable X^ after n trials. 
If the state at a trial n is i, the probability that the next state is 
equal to j is defined by Eiiia 

-P(^n+l ~ j\^n — '^ni ^n—1 ~ '^n—1') • • • : ^0 ~ '^o) 

= P{Xn+l = j\Xn =i) = Pij 

Equation ([^ is known as the Markov property and says that 
the state of a random variable after n + 1 transitions, only 
depends on the state of the random variable after n transitions. 
In other words, the future of the system will only depend on 
the present state ElISl. 

If we want to know the probabilities of a single state it is 
useful to use the total probability theorem which is given by 


ca 

P{B) = J2PiB\Ai)P{Ai) (8) 

i 

Where the index i runs for all the states Ai that B depends 
on. If B is the event of = j then P(Xn+i = j) = 

Pj{n + 1). Ai is the event of X^ = i, thus P{Xn = i) = 
Pi (n), then the theorem transforms into the following equation 

il[l3 

Pj{n + 1) = ^pij{n)pi{n) (9) 

i 

This notation allows us to interpret equation ([^ as a multi¬ 
plication of a horizontal vector formed with the probabilities 
of some certain state with a matrix formed by the pij . So if 
there are N possible states, p(n) = [pioin), ... 
and W is the matrix with elements pij , the relation between 
states can be written as CSIEI 

p(n + 1) = p(n)W(n) (10) 

Given the initial condition p(0) we can determine the solu¬ 
tion to Eq. 

III. THE MODEL 

Suppose that we are given a collection of N infinite charged 
planes, half of them have a constant charge distribution a and 
the other half have a constant charge distribution —a. We can 
write this condition as 

N-l 

J2<^i^n)=0 (11) 

n=0 

If we randomly place the N infinite charged planes par¬ 
allel to each other along the 2 : axis at position Zn = n e 
{0,1,..., — 1}, there is no way to know the values of 

the electric field in between the planes unless we measure the 
charge on the planes. Instead of that we will use Markov chain 
theory to analyze the problem. 

We will first consider a system of four charged planes. 
E{z) equals zero at the left and right side of the configuration 
due to the neutrality of the system HI . After crossing the first 
charged plane, E{z) would increase or decrease in an amount 
of cr/eo depending if the first charged plane had a positive 
or negative charge, and the same for the remaining planes 
according to equation ([^. In figure there is an example of 
how the electric field behaves for a known charge density. 
In figure we show the possible values the electric field can 
take when the charge density is not known. 

Eor the sake of simplicity, we normalized the electric field 
dividing by cr/eo so 

E{z) = 


( 12 ) 
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Due to this normalization E{z) takes only integer values i. 
Defining Ei (n) as our random variable, where n is the region 
between Zn and Zn+i, our state space S is defined by integer 
values between the interval [—N/2^N/2]. To calculate the 
transition probabilities pij we use the fact that the boundary 
condition in equation ([^ prevent us to have a value bigger 
than i + 1 or smaller than i — 1 at a place z = n + 1 when 
Ei{n) = i. This results in pij = 0 for j > i + 1 and for 
j < i + 1. After crossing a charged plane the field has to 
have a different value, so pa = 0. The last values for pij are 
determined by the probability of finding a negative or positive 
charged plane at the transition point. Defining W± (n) as the 
probability to find a positive or negative charged plane at a 
place Zn = n, we get the following results for pij 


f W-^{n) if j = i + 1 

Pij = I W-(n) (13) 

[ 0 otherwise 

A simple way to calculate the factors W± (n) its by dividing 
the remaining positive or negative planes by the remaining to¬ 
tal planes. Let p± be the number of positive (negative) planes 
between 0 < z < n, then = n. Also the value of the 

electric field i at a point n is completely defined by the amount 
of positive and negative plates that have left behind, therefore 
= i. Using these two relations we can solve for 


P+ = 


n i 
2 


(14) 


The remaining positive planes are N/2—p-^ and the remaining 
total planes are N — n, then 1U+ is given by 


W+(n) = 


N _ 2+n 

2 2 


N — n 



i 

N — n 


) 


(15) 


andiy_(n) = l-W+{n). 

We known that the system is neutral, so the electric field 
must be zero at 2 : = 0. In other words the probability for 
E to be zero at z = 0 equals one, and the other cases have 
0 probability. Using our notation the last statement can be 
written as po{0) = 1 andpi(O) = 0 for i ^ j. 

With all these information we can use Eq. (1^) to obtain 
the most probable value of the electrostatic field in between 
the charged planes. 
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IV. RESULTS 

In this section we present our results for a the system which 
consists of 100 charged planes (N = 100).Figure|^ shows a 
3D plot for p(n). Note that the highest values for the prob¬ 
ability are reached at n = 0 and n = N — 1 with i = 0, 
which shows the fact that the electric field must be zero at 
these points. In the direction of the increasing n the proba¬ 
bilities spreads out on the possible states, so the electric field 
becomes more uncertain at places near the middle of the con¬ 
figuration. After that point, the probabilities change to finally 
reach 1 at n = — 1. 

Its important to note that these probability functions for 
each transition are discontinuous functions. The plot in figure 
[^represents a cross section for p(n) at the middle of the 
configuration (n = 50). In this plot it is easy to see that there 
are jumps between probability zero and nonzero probability 
along the possible states. This is due to the boundary 
condition given in Eq. 0- Note that the field cannot remain 
the same so there are forbidden states at a certain transition n. 

Now we can calculate the two physical quantities that we 
are interested in which are the mean electrostatic force and 
mean electrostatic energy density. It is well known that the 
electrostatic force is given by F = qE and the energy density 
hy 6 e = eoE‘^{z)l2. In order to calculate these two quanti¬ 
ties we need to calculate the first and second moments for the 
electrostatic field, i.e. {E{n)) and (F^(n)). To calculate such 
moments we use the ^-th moment formula EllISl 


(E'=(n)) = ^z''Pi(n) (16) 

i 

Where the index i runs through all S. In order to calcu¬ 
late these moments let us rewrite the form of equation ( p^ as 
follows. Most of the probabilities pij equals to zero, so the 
probability pi(n + 1) can be calculated adding the 2 possible 
paths from the last transition given by the possible values at 
the transition n. Each of the paths may me multiplied by the 
factors W±{n) |[T3]| as is depicted in figure So the new 



FIG. 2: Possible values for E(z) inside a four charged planes system 
FIG. 1: Values for E(z) inside a given configuration of charged planes where the charge distribution is not known 
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FIG. 3: State probabilities for every transition 



FIG. 4: State probabilities after 50 transitions 


obtain 


N/2 

i=-Nl2 

N/2 

i=-N/2 


(18) 


Its important to note that PAr/ 2 +i(^) = P-Ar/ 2 -i(^) = 0 
because N/2 + 1 and —N/2 — 1 does not exist in S. Also 
W± {±N/2^ n) = 0 because the field cannot get a value higher 
(lower) than ±N/2. Using these facts and setting k = 1 in 
GU we obtain 

{Ein + l)) = ^~//\ E{n)) (19) 


Due to the neutrality of the system we can write (^(0)) = 0. 
Using this initial value and substituting it in ( p^ its obvious 
that {E{n)) = 0 for all n inside the plane configuration as 
expected ca. Therefore the mean force is zero everywhere. 

Setting k = 2m ( p^ and doing a little bit of algebra, equa¬ 
tion turns into 

_ 7\r _ _ p _ 

+ 1)) = N-n ^ 


Let A{n) be the factor that multiplies {E‘^{n)) in equation 
(19), and knowing that (^^(0)) = 0 and using recursively 
the difference equation in ( [20| ) ifT^ . the states for the second 
moment are given by 


form of equation ( p^ is given by 


Pi(n-hl) = Pi+i(n)VF_(i-h l,n) 
+ Pi-i(n)VL+(i - l,n) 


(17) 


(^ 2 ( 1 )) = 1 

{E\2))=A{l) + l 

{E^{2,)) = A{2)A{l)+A{2) + l 

( 21 ) 

{E^{n))=lN^ilA{n-j) 


Multiplying both sides of equation ( [T^ by and summing 
from —N/2 to N/2 we may use the definition given in ( p^ to 



Substituting A{n) in equation (21), the second moment is 
given by 


n—1 


/p2/ 1 I {N-n){N-n-l) 

^Ei^N-n + i){N-n + i-l) ^ ^ 

To see how the energy evolves along the z axis we used the 
same example with = 100 as before and numerically calcu¬ 
lated the energy using ([22|).The result is depicted in figure 
The the mean energy may be calculated using the definition of 


our random variable given in (12) as 


{5^) = ^{E\n)) =—{E\n)) (23) 


FIG. 5: Relation between transition n and n + 1 


We can use standard probability theory to validate our calcu¬ 
lations for the example given above. For the case of = 100 
there is a probability of 50!50!/100! that the electrostatic field 
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FIG. 6: Evolution of the second moment along the z axis 


reaches the maximum (minimum) value ±50 after 50 transi¬ 
tions. Using our approach this means that the probability to 
get a value of 50 is obtained by 


50 


P5o( 50) = IIW"+(50-i,50-i) 


(24) 


i=i 


Equation ( [24| ) means that we need 50 positive charged planes 
in order to reach the maximum value of the electrostatic field. 
Using the fact that lU+(i,n) = {1 — i/{N — n))/2 we get 
lU+(50 — j, 50 — j) = j !(50 ± j), inserting this into Eq.(24) 
we get 


„ fr ^ (l)(2)(3)---(49)(50) 50150! 

Psolou; Ugo + j (51)(52)(53)---(99)(100) (100)! 

(25) 

Equation ( [25] ) shows that our result is consistent with the stan¬ 
dard probability theory. 


V. CONCLUSIONS 

We have shown that it is possible to perform a random walk 
with and electrostatic field by means of several parallel infinite 
charged planes in which the surface charge distribution is not 
explicitly known. We have worked out the special case where 
the charged planes have a constant surface charge distribution 
±cr and the overall electrostatic system is neutral, i.e. there is 
the same number of positive and negative charged planes. We 
use Markov chain theory in order to give the most probable 
value for the electrostatic field in between the charged planes 
and use these results to obtain the mean electrostatic force and 


energy density. Our model gives rise to a stochastic law when 
the charge distribution is not deterministic. 
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